Heat-kernel calculation of quark determinant 
and computer algebra 



A.A.Bel'kov^ D.Ebert^ A.V.Lanyov^ A.Schaale^ 

^ Particle Physics Laboratory, Joint Institute for Nuclear Research, 
Head Post Office, P.O.Box 79, 101000 Moscow, Russia 
^ Institut fiir Elementarteilchenphysik, Humboldt-Universitat, 

Invalidenstrafie 110, O-1040 Berlin, Germany 
^ DESY-Zeuthen Platanenallee 6, 0-1615 Zeuthen, Germany 

February 1, 2008 

Abstract 

In this paper there we describe the calculational background of deriving a strong meson 
Lagrangian from the Nambu-Jona-Lasinio quark model using the computer algebra systems 
FORM and REDUCE in recursive algorithms, based on the heat-kernel method for the calculation 
of the quark determinant. 







Computer Algebra Systems (CAS) such as FORM [Q or REDUCE have been successfully used 
in the high energy physics for a long time, especially in the field of the perturbative theory calculations 
of higher order contributions in the standard model. Besides these "classical" fields of application of 
CAS, they can be also used for solving various problems connected with non-renormalizable models 
such as chiral meson models. In the following we describe the calculational background of the recent 
derivation |^ of meson Lagrangian describing strong interactions of scalar, pseudoscalar, vector and 
axial- vector fields. 

The strong interactions of quarks at high energies (i? ^ 1 GeV) are described by the Quantum 
Chromodynamics (QCD) Q. The QCD Lagrangian has following form: 

with 

Gf^u — dfj^Gu — diiG^ + ig\Gfi, Gy\, 
D = ^^{d,-igG,), G, = Y^Gl^, 

a 

where q is the quark field, G^ is the giuon field, g is the coupling constant, mP^ is the current quark 
mass, are the Dirac matrices, Ac the SU (3)c matrices. 

Some special properties of the QCD like the confinement of quarks and the momentum dependence 
of the strong coupling constant as(9^) which becomes large at low momenta {a^ « 1) make it 
impossible to use the QCD for an adequate description of meson processes at low energies. 

Until now it was not possible to derive a meson theory from the QCD. On the other side, meson 
processes have been studied in a framework of chiral Lagrangians since 1960s. One possible way to 
built up a QCD-based chiral meson theory is to start with some effective quark Lagrangian, e.g. with 
the Nambu-Jona-Lasinio model (NJL) 

The NJL Lagrangian of the effective four-quark interaction has the form: 

-Cnjl = q{id - mo)q + dnt (1) 

with 

A„t - 2Gi ^{q\\)q? + {qil'\\)qf'^ - 2G2 l^{qr\\)qf + (97 V ^A}?)^ | , 

Gi and G2 are some empirical constants. 

Let us show that the NJL model can be motivated from the nonperturbative QCD using some 
assumption about the infrared behavior of gluons. The generating functional of QCD has the form 
i: 

qi{iB^m))qf +gqfG''^qf - -^G^'G^, 



J VqVqVGexpi^i J d^x 
After the integration over the gluon fields one will get 



Z./l„P,-exp{./...«.S-.»?),}.xp(.1'M) 
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with 

as flavour-singlet quark current and W becomes 

71=2 ■ i=l 

The functions i*^! '^" contain all information about the gluons. An analytical integration is not 
possible because the behavior of the gluon propagator and the running coupling constant at long 
distances are unknown. 

Neglecting higher order Green functions one can get some local approximation for the gluon 
propagator D^^{x, y) oc g^,y5^{x — y) and find some theoretical base for the approximation .Cqcd ^ 
£njl- In this approximation it is impossible to predict the values of the constants Gi and G2 from 
QCD principles. They are to be determined experimentally. 

One can see that the group structure of the QCD 

5[/(3)J5'=^' ® SU{Nf)L ® SU{Nf)R ® U{1) ® S 
and that of the NJL model 

^j^^g^giobai ^ sU{Nf)L ® SU{Nf)R ® U{1) S 

are very similar {S presents the set of discrete symmetries as C, P and T conjugation). 

After this introduction of the NJL model we will now derive a meson Lagrangian from the NJL 
Lagrangian (|l|) using path integral techniques. First one has to introduce collective fields with the 
form: 

which will correspond later to scalar (S), pseudoscalar (P), vector (V) and axial (A) colourless fields. 
The following notations have been used: 



Now the Lagrangian can be redefined: 

^NJL = ti''^^'^ " 4^ *^ ^^i^ + +q{id + V + - Pii{(j) + mo) + Pl(<A+ + mo)} q 



with (j) — S + iP, V — V^^^, A — A^7,i and Pr/l = 5(1 ^ 7^) ^ chiral projectors. This substitution 
transforms the Lagrangian into a bilinear form in quark fields. Now the integration over quark 
fields becomes possible and the generating functional of the NJL model is 



Z = j Vc^Vcji+VVVA exp ijd^x-^^ tr 0+0 - ^ tr [V^ + 4)) 
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with 

Zf{(l),<j)+,V,A) = J VqVqexp{i J d'^xqiBq) = (detiD) (2) 

and 

if> = id + V + Aj'^ - PbU + mo) + PL{(j)'^ + mo) 

= i{d + Ar) - (0 + mo)PB. + i{d + Al) - (0+ + mo)PL 

as Dirac operator. Ar/i^ = V ± A are right- and left-handed vector fields. The quark determinant 
(|^), describing the interactions of mesons, can be evaluated in quark loops |Q, ^. 

At this place it makes sense to refer to some different ways in the calculation of the quark- 
determinant (H). In contrary to the "straight" method of calculating the quark loops [Q we have 
used the heat-kernel method . Of course, the physical results must not depend on the method of 
calculation but there are some technical reasons for using this method here. The main advantage of 
this method is that the recursive algorithms of this method can be adapted on CAS quite effectively. 
The infinite contributions can be conveniently separated too. 

The absolute value of the quark determinant (H) in "proper-time" regularization is defined as 

log|detiD| = --Tr' log(D+D) = -- / dr- Tr' exp(-D+DT) (3) 
2 2 7i/A2 T 

with A as regularization parameter. The trace Tr' is to be understood as a space-time integration 
and a "normal" trace over Dirac, colour and flavour indices: 



Tr' = j d^xtr', tr' = trc • try • tr-^ 



The operator D+D can be presented as 

D+D ^ d^df" + a{x) + /x^ 

with 

dp = 5p + rp, T^ = V^, + A^^^, a{x)^^VH + H+H + ^[J^',-f'']T^,~^J.^ 

Here fi plays the role of some free parameter which will fix the regularization in the region of low 
momenta. Below it will be identified to the constituent quark mass. We used the following notations: 

H = Pii{<i) + mo) + Pl(0+ +mo) = S + i-f^P, 

r^. = K, d.] = a^r, - d,T^ + [r^, r,] = f^^ + j'f;^,, 

with i^^^ as field strength tensors 

F^u = d^A, - d,A^ + [V^,A,] + [A^, K] 



3 



and 



V^i/ = a^i? + [V^, H] - 7^{A^, H} 



as covariant derivative. 

The main idea of the heat-kernel method is to evaluate 



< X I exp(-D+Dr) | y > 



around its nonperturbated part 



< X I exp(-(n + /i^)r) I y >= 



(47rT)2 



in powers of proper-time r with the so-called Seeley-deWitt coefRcients hk{x, y) 



< X I exp(-D+Dr) | y > = 



(47rT)2 



-M^r + (K-y)V(4 



After integration over r in one gets the following expression for log | det iT) 



with 



^ I„g(*. D-D). 4 Jip 



r(Q;, x) 



t .Q-l 



2k 



Tr 



as incomplete gamma function. Using the definition of the gamma function one can separate the 
divergent and finite parts: 

1 



- log(det D+D) = Bpoi + Slog + %n- 



Here 



R =1^ 
Pol 2 (4^)2 



4 -[^ 4 



has a pole at a; = IJ- /A — 



1 1 



rr(0,a;) 



log 2 (4^)2 
is logarithmic divergent. The finite part has the form: 



i/i'^ Tr' ho - Tr' /ii -I- Tr' /12 



1 1 



fin 



2 (47r)= 



^ /i4-2'=r(A: - 2, x) Tr' /i^ 



fc=2 



The main technical problem now is to calculate the Seeley-deWitt coefficients hk- One can find the 
calculated heat coefficients up to fc = 3 |^ and with simplifications up to fc = 6 in |^ . In this paper we 
have calculated the full coefRcients up to the order n = 4 and and also present the minimal parts of 
the heat-coefficients h^^. Let us show the general heat-kernel techniques in detail: The determinant 
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of the positive definite operator A is defined in proper-time regularization by the following integral 
relation 



log(detA) = - / —Tr:'K{T), 



1/A2 

where K(t) = e~^'^ is the so-called "heat kernel" which satisfies the equation 

^K{t)+AK{t)=Q 

with the boundary condition 

K{t = 0) = 1. 

In the case discussed in the paper the operator A has the structure 

A = d^d^" + a{x)+ 

where 

contains an differential operator and a{x) docs not. 

The asymptotic behavior of A at long distances is corresponding with the infinite part of log(det A) , 
which is defined by the "free" part 

Ao=n + /x2, nsS^a". 

Using the ansatz 

K = KoH 

it is convenient to separate the "free" part Kq from the heat kernel. In coordinate representation Kq 
reads 

ifo(x,y;T) =<x\K^\y>= ^_L_e-'''-+(--'')'/(^-) 
and also satisfies the equation 

^+AoKo = 0, 

with the boundary condition 

Kq{t = 0) = 1. 

The "interaction" part H of the heat kernel satisfies the equation 

d 1 

iir + -^t^d^ + d^d^ + a)H{x, y; r) = 0, 

OT T 

H{x,y = x:T = ()) = l, (4) 

where = x^ — y^ and the differential operator d^ acting only on x. Using now an expansion for 
H{t) in powers of r 

oo 

fc=l 



H{x, y-r) =^ hk{x, y) ■ 
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one will get from (0) the recursive relation 



(n + 1 + z,,d^)/i„+i + (a + d^,d^')K = 

with boundary condition 

z^#/io = 0. (5) 

The heat coefficients hk{x) = h{x,y = x) are defined by the diagonal part of the recursive relation 
i) 

nhn{x) = -[d^d^/i„_i(a;,y)] \y=^ -ahn-i{x) (6) 

which gives for n = 

hi{x^ x) = —{d^df'^ + a)ho{x, x). 
Here we use the more convenient notations 

da... dphn{x, x) = [da... dphn{x, y)] \y=^ . 

The expression can be trivially calculated once the {n — 1)*'* order heat coefficient is known. To 
find the derivative terms d^d^hn-i one applies the operator d^d^^ on the recursion relation (^) for 
n — \ and puts afterwards y = x. This introduces derivatives up to the order of hn-2, and in addition 
derivatives of a. One continues to apply repeatedly operators d^i on the recursion relation (^) for 
smaller and smaller n until the desired n*'* order heat coefficient is completely expressed in terms 
differing only by a permutation of the cZ^'s. To get an explicit expression for the multiple derivatives 
of hn one has to bring the d^'s in all of these terms into the same order. The reordering of the d^'s 
introduces multiple commutators of the form: 

Kfj,i, = [dfj,,d^] = Tfj,^, K\f_,y = [dx^K^i,], K^x^^ — [dnjKx^ij] , etc. (7) 

Similarly one has to bring the a's to the left of all d^t's. This introduces the multiple commutators 

Sf_i^[d^,a], Sf_i^ ^ [d^,Si,], Sxtj.iy ^ [dx, Sf_ii,] , etc. (8) 

In this way the heat coefficients can be completely expressed in terms of multiple derivatives of Hq, 
of multiple commutators (|^) and (||), of lower-order heat coefficients, and of a's, with the latter ones 
being always on the left of the derivative operators d^. The multiple derivatives 

da. . . d^ho{x) = da - . . d^ho{x, y) \y=x 

are finally calculated by applying the d^'s repeatedly on the differential equation (||) defining /iq. 

Following the strategy outlined above the calculation of the heat coefficients is straightforward 
but cumbersome. The very lengthy calculations can be performed only by computer support. The 
calculation of the heat-coefficients is a recursive process which can be done by CAS very conveniently. 

For the calculations of the heat coefficients we used the CAS REDUCE and FORM. The reason for 
the use of different CAS for the solution of one problem is connected with the special features of the 
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systems. FORM is designed for doing fast arithmetic with large formulae and a special attention was 
to implement there structures commonly used in high energy physics, e.g. summation over dummy 
indices. It easy easy to control even difficult calculation processes step by step since operations are 
only done on request. However, the scope of FORM operations is limited and if some particular sphere 
of symbolic calculations is somewhat biased from the generally used area of multiloop calculations or 
Dirac algebra, then a need for unimplemented operations can arise. In our case of calculations with 
nonrenormalisable chiral models in low-energy meson physics there is a need not only for summation 
over some dummy indices but also for the following operations. 

First of all one needs the operation of cyclic shifting of non-commuting operators products under 
the trace sign. This operator is required since the following parts of Lagrangian are equivalent: 

£1 = tr (9^$ • $ 9^$ • $) and £2 ^ tr ($ 5^$ ■ $ 9^$) 

and one should transform them to some unique form, either the former or the latter one. 

Secondly, one way of testing the obtained physical results is to test their symmetry properties. In 
our case we have Hermitian conjugation and the transformation of intrinsic parity. One needs some 
means to test automatically the symmetry properties of the results. 

Both FORM and REDUCE lack the described operations. FORM has no possibility of adding new 
non-trivial operations. For this purpose we used CAS REDUCE which is a completely open system 
allowing the user to access directly the internal data structure and implement new operations and 
data types by the "symbolic" style programming. We implemented the required operations absent in 
FORM and also such as summation over dummy indices. In such a way the whole work described in 
this paper could be done only with the help of the single REDUCE system. Let us briefly describe 
the corresponding REDUCE functions from our package |^ for applications in meson chiral models. 
Each REDUCE example below follows the usual mathematical notation. 

traceshif t (X) - transforms expression to the unique form using the possibility to perform cyclic 
permutation of the products of matrix symbols under the trace operation, e.g. 

tr (9^$ • $ 9^$ • $) — > tr (9^$ • $ 9^$ • $) 

traceshif t(fi(mu)*fi()*fi(niu)*fi()) > f i (mu) *f i () *f i (mu) *f i () 

tr ($ 9^$ • $ 9^$) — > tr (9^$ • $ 9^$ • $) 

traceshift(fi()*f i(mu)*fi()*fi(mu)) > f i (mu) *f i () *f i (mu) *f i () 

order ind(X) - transforms expression to the unique form using the possibility to redesignate the 
dummy indices, e.g. 

tr (9^$ • $ 9^$ • — > tr (9^$ • $ 9^$ • 
orderindCf i(mu)*f i()*f i(mu)*f iO) — > f i (mu) *f i () *f i (mu) *f i () 

tr (9^$ • $ 9^$ ■ — > tr (9^$ • $ 9^$ • 
orderindCf i(nu)*f i()*f i(nu)*f iO) — > f i (mu) *f i () *f i (mu) *f i (). 
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After voluminous computations one gets the complex expressions for heat-coefficients hi,. . .hi: 

hoix) = 1, 
hi{x) = —a, 

tv'[h,{x)] = tr'|^(r^.)2 + ia2|, 

tr' [hs{x)] = -Itr' |2a3 - S^S'^ + a{r,^)' - ^{K^p.f - ^(i^^a/s)'} , 
tr' [h,ix)] = tr' I la^ + ^a^^^ + aS.Si^ + ^7(3%)' - {S^^f 

(Qo ^ 1 27 1 \ 
+ K^'f^i^'' K°""''^ K^'"'^" 
25200 1575 5040 600 / 

13 47 17 

+ 12600^^''^''"^''"^"'' 16800 + 25200 

+ T575^-'^'*"-'^""^' 25200 ^""'""^"'^ ' ~ 12600 ^^"'^""^^ ^ 1575^^''"""^^ 

6300 ^ 5600^ ' 5040 



1800 ^ 25200 

In the same way the next orders of heat expansion coefficients hi can be obtained using the 
developed calculation technique based on the usage of computer algebra. For simplicity wc present 
below expressions only for minimal parts of heat-cocfficicnts, i.e. only for the parts which do not 
vanish in the pseudoscalar region of the theory when — = 0; 

tv'[h,{xrn = -tr'l^^a\a' + S,Sn + ^a'S,^ + 2{aS,f 

+ ^ [I0a5^(5^^ + S,--^) - 2a{S^,f + na{Sp' + S>'\S^ + SaS^'^'^S,] 
11 1 Q 2 



1008 ^ 2800 225 

1 

25200 
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( 23 937 23 \ 1 

4800^"^'' + 302400'^''^'''^ + 10800 ) + 252 



352800 A^-v « ■ ■ " V3360 ^ 11340 

17 1 
^226800''^^''""^' ^ 317520" + S^''^^"'' + 775''''^"a) 5, 

~30^ [53(5^5")^ + {S.S.f] + ^^S,S. (1575''„«^ + 2985''''«„) 
+ 70^0^"^" (315„«^'' + 585^''«„ + 475^«„^) + 
+ 14112^''^'^^""''' + 10^^'' ^^""^ ^^^^""^ + ^^^-""^ + 355-^5,"J 
+ 21168^'* I^''''^""- + S^^i5S•'^^" + 25^"")] + {S^^^ + S,-"^) 



/ 83 „„„ 1 „ 1 



Vl41120'^'"'" + 9408^"" J + 30240'^'' + ^^'^""'^^ ^'''^ 

^- — (701(5„'')^ + 5835''''5„„S'?) - S^'JS^^f 

635040 ^ V M ; 316386 ^ ^ 

2 „„ 1 



Son/ acoc/^n (619*5' uS^a ~i~ ISOtS^a^* jy) 



2835-^ ^ 952560 ' ' 

226800^ '^""''^ 12700800^ " ^ 66150 ^ 
1 13 1 

c/i 1/ qa)} qui^a q -0 ( 



52920 604800 

To obtain these expressions for the heat coefficients we have extensively made use of the cychc 
properties of the trace and of the Jacobi identity. To obtain effective meson Lagrangians in terms of 
collective fields one should calculate in tr' hi{x) the trace over Dirac indices. 

We would like to thank V.P.Gcrdt for useful discussions. Two of the us (A.A.Bel'kov and A.V.Lanyov) 
are grateful for the hospitality extended to them at DESY-Zeuthen. 
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